Fronts with a Growth Cutoff but Speed Higher than v* 
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Fronts, propagating into an unstable state = 0, whose asymptotic speed u as is equal to the 
linear spreading speed v* of infinitesimal perturbations about that state (so-called pulled fronts) are 
very sensitive to changes in the growth rate f(4>) for (f> <C 1. It was recently found that with a small 
cutoff, f(cf>) = for (j) < e, « as converges to v* very slowly from below, as In -2 e. Here we show that 
with such a cutoff and a small enhancement of the growth rate for small <j> behind it, one can have 
n as > v* , even in the limit e — » 0. The effect is confirmed in a stochastic lattice model simulation 
where the growth rules for a few particles per site are accordingly modified. 

PACS numbers: 05.45.-a, 05.70.Ln, 47.20.Ky 
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Pulled fronts are those fronts that propagate into a lin- 
early unstable state, and whose asymptotic front speed 
w as equals the linear spreading speed v* of infinitesimal 
perturbations about the unstable state ji], ^ The 
name pulled front refers to the picture that in the lead- 
ing edge of these fronts, the perturbation about the un- 
stable state grows and spreads with speed «*, while the 
rest of the front gets "pulled along" by the leading edge. 
That this notion is not merely an intuitive picture but 
can be turned into a mathematically precise analysis is 
illustrated by the recent derivation of exact results for 
the general power law convergence of the front speed 
to the asymptotic value v* j|. Fronts which propagate 
into a linearly unstable state and whose asymptotic speed 
^as > v * are refered to as pushed, as it is the nonlinear 
growth in the region behind the leading edge that pushes 
their front speed to higher values. If the state is not lin- 
early unstable, then v* is trivially zero; in such cases the 
front propagation is always dominated by the nonlinear 
growth in the front region itself, and hence fronts in this 
case are in a sense "pushed" too. 

For the field 4>(x,t), the dynamics of fronts that we 
consider in this paper is given by the usual nonlinear 
diffusion equation 



d(f> d 2 (j) 



(i) 



In the standard case, the growth function f(6) has the 
form f{4>) = (j) — <j) n , with n > 1. Equation ([!]) has two 
stationary states for d>(x,t): <j>(x,t) = and (f>(x,t) = 1. 
Of these, <f>{x,t) = 1 is stable and <f>(x,t) — is unstable. 
The asymptotic speed of (pulled) fronts propagating from 
(f>(x, t) = 1 into <f»(x, t) = in Eq. (g) is v* = 2. 

The sensitivity of pulled fronts to the precise dynam- 
ics for small perturbations about the unstable state has 
recently surfaced in a remarkable way |jj . Often, in equa- 
tions like (g), the field <f)(x,t) is the density of particles 
in a continuum description. If one then considers fronts 
in stochastic particle model versions of (g), the linear 
growth term in f(4>) implies that for small particle den- 
sity, the rate at which new particles are created is propor- 
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FIG. 1: (a) Shape of the function f(cj>, e) used by Brunet 
and Derrida to study the effect of a finite particle cutoff in 
the growth rate on the front speed, (b) The growth function 
f(4>, e, r) (thick line) we analyze in this paper. In both cases 
we have kept only the linear term of /(</>) to plot the graphs, 
since e <C 1, so that the nonlinear terms in / are much smaller 
than the linear terms. 



tional to the density itself. Brunet and Derrida || were 
the first to realize the fact that for new particles to be 
created in any given realization, the density must be at 
least one "quantum" of particle density strong, and that 
this provides a natural lower cutoff for the growth that 
strongly affects the front speed. Indeed, to mimic this 
effect, they considered a deterministic front of the type 
in Eq. ([l]) with n = 3, and by hand introduced a cutoff 
of the type sketched in Fig. ga in the growth function at 
cj) = e <C 1. In this paper, we denote their growth func- 
tion by f(<j), e) = [(f) — 3 ] 0(0 — e) where O is the unit 
step function. For small e, the asymptotic front speed 
w as (e) was then found to be [|| 



Was(e) 



(lne) : 



+ 



(2) 



Brunet and Derrida subsequently identified e with 1/N, 
where N is the average number of particles at the sat- 
uration state of the front, corresponding to the stable 
state <j)(x, t) = 1 of the density field. The slow loga- 
rithmic convergence to the asymptotic front speed from 
below as a function of N, implied by Eq. (0), has been 



2 



confirmed in various studies of stochastic lattice models 
d 1, |, 0|, §• Note that for < e, the growth function 
/ vanishes, and as a result, strictly speaking, the state 
<f> = is not linearly unstable; hence fronts in this model 
are always weakly pushed for any nonzero value of e |l0| . 

In this paper, we demonstrate an even more surprising 
aspect of the sensitivity to small changes in the growth 
function / of the "pulled" fronts, that we have at e = 0: 
if / is sufficiently enhanced in a range of <f> of the order 
of e, the asymptotic front speed v as can become larger 
than v* and not converge to v* as e — > 0. For fluctu- 
ating fronts, this implies that if the stochastic growth 
rates for small occupation densities rii are somewhwat 



enhanced over a linear behavior 



then such stochas- 



tic fronts may move faster than v* and never converge 
to their naive mean field limit for N — > oo. This effect 
may be of relevance for the coarse-grained field theory for 
DLA, as it is empirically known to be essential to modify 
the growth function for small cluster densities [H . 

We now discuss our results first, and then summarize 
their derivation. 

To be specific, we consider the nonlinear diffusion 
equation (Q) with the growth function sketched in Fig. |l|b, 



f((j),e,r) = /(</>,e) for <f) < £ and for 
= e/r for e < cf) < e/r , 



> e/r 



(3) 



lim f(<f>, e,r) = f{4>) , 

e— >0 



with r < 1. We show that while for any fixed value of r, 

(4) 

the asymptotic front speed v as (e, r) has the property that 

for r > r c , 



where 



lim u as (e, r) = v* 
lim v as (e, r) > v* 

e— »0 



■ («*»- a) 



for r < r c 



0.283833. 



(5) 



(6) 



Hereafter, to save writing, we denote u as (e, r) simply by 
v. For e — * 0, the asymptotic speed at a given value of 
r < r c in our model is given by the relation 
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from which the value of r c , given by Eq. (g), follows. 

These expressions show that the limits do not commute 
for r < r c : taking the limit e —> first in / yields a front 
speed v* but the limit v as (e — » 0, r) > w*. The reason 
is that for r < r c there is always a little tail of the front 
that runs faster than v* and makes 6 nonzero. Once 6 is 
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FIG. 2: Comparison of simulation data for v as (e,r) with the 
analytical prediction (Q), which is plotted as the solid line. 
The solid dots in Fig. ti represent the numerical data for 



Eq. 



with e = 2 x 10 and n = 3. The crosses are the 



data points for fronts in the stochastic growth model described 
in the text. Inset: illustration of the leading order rate of 
convergence of the u as (e, r) curve to the e — » limit, by means 
of the schematic dashed curve. 
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FIG. 3: Numerical data for v as (e, r) as a function of e for n = 
2 and 3, at r = 0.2. The graph demonstrates the insensitivity 
of v to e for small values of e (note the fine scale on the vertical 
axis), as well as the convergence as e n_1 to its e —* value 
~ 2.246, given by Eq. (^), for two different values of n. 



nonzero, growth continues and the region behind it just 
has to follow it with the same asymptotic speed. 

Our analysis is corroborated by numerical results ob- 
tained by solving Eq. (|l]) forward in time, (with Gaussian 
initial conditions). The data for v vs. r at e = 2 x 1CP 5 
are shown as solid dots in Fig. 0. Note that for r < r c , 
the solid dots fall on top of our prediction drawn 
with a solid line, while for r > r c , they systematically 
fall below the solid line v = v* . The reason for it is the 
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difference between the rates of convergence as e — > 0, 
which is illustrated in the inset of Fig. || by means of the 
schematically drawn dashed line. The arrows in the inset 
indicate the rate of convergence of the dashed r-v curve 
towards the limiting one, given by Eq. (Q). For r > r c , 
the convergence is ~ In -2 e like in the case for r = 1, 
analyzed in Ref. ||]; but for r < r c the convergence is 
much faster, e" _1 . This latter behavior is illustrated for 
r = 0.2 and n = 2, 3 in Fig. || — note the fine scale on 
the vertical axis! 

That the effect of increasing asymptotic speed with de- 
creasing r below r c is a real effect for stochastic fronts too 
is illustrated by the crosses in Fig. ||: these represent the 
data for the average speed of fronts in a reaction-diffusion 
system 2X, for discrete X particles on a lattice with 
N = 10 4 Jl^], where the growth rates have been modified 
when the number of particles rij on a lattice site i is less 
than 1/r. In accord with the shape of the growth func- 
tion / illustrated in Fig. [j]b, the rate at which particles 
are created at a lattice site i with 1 < rij < 1/r particles 
is simply taken to be the same as the rate for rii = 1/r 
(corresponding to the integral values 1/r = 1, 2, 3, 4, 5, 7 
and 10, due to the discreteness of particles). As one can 
see from Fig. already when r — 0.5, i.e., when only the 
growth rate at lattice sites with one particle is increased 
by a factor 2, the asymptotic growth speed is above the 
value v* = 2. 

In the remainder of this paper, we derive the analyti- 
cal results for the nonlinear diffusion equation with the 
growth function (|3|). Our analysis is based on the fol- 
lowing observation: for e = 0, it is well known that the 
nonlinear diffusion equation allows a continuous family 
of front solutions with v > v*. When such fronts so- 
lutions are parametrized by their velocity v, and when 
the growth rate is modified to allow a transition to a 
"pushed" front with velocity , it is also known g, |j| 
that solutions with v < are unstable to a localized 
mode. In our analysis, we therefore consider a front with 
a given fixed velocity v and, for small e, determine when 
upon decreasing r a localized mode of the stability oper- 
ator crosses the eigenvalue zero. In the limit e — > this 
marks the selected pushed front in the r-v diagram. 

To carry out the linear stability analysis of the front 
solution, it is convenient to follow the standard route 
of transforming the linear eigenvalue equation into a 
Schrodinger eigenvalue problem |l|, We consider a 
function 4>(x, t), which is infinitesimally different from the 
asymptotic front solution </> as (C) in the comoving frame 
£ = x — vt, i.e., (f>(x,t) — </> as (£) +v(^t)- Upon lineariz- 
ing Eq. (Q) in the comoving frame, one finds that the 
function 7/(£, t) obeys the following equation: 



dr/ 



dr/ 



(8) 



Since this equation is linear in 77, the question of sta- 
bility can be answered by studying the spectrum of the 
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FIG. 4: The potential V(£) for v > v* and infinitesimally 
small e in the Schrodinger operator that determines the tem- 
poral eigenvalues of the stability analysis. £2 marks the posi- 
tion of the region of finite width where the potential crosses 
over from the asymptotic value on the left where (j) as w 1 to 
the value in the well where <^) as <C 1, £1 the position of the step 
and £0 the position of the delta function term in the potential. 



temporal eigenvalues. To this end, we express ?y(£, t) as 



= P - Et f>-</ 2 



(9) 



which converts Eq. (||) to a one-dimensional Schrodinger 
equation for a particle in a potential with h 2 /2m = 1: 



d£ 2 + 4 



MO = EMO, (10) 



In Eq. @, the quantity V(£) = 



plays the role of the potential, 
plicitly from the expression (0) 



v 

T 



It is easily obtained ex- 
for f{4>, e, r) as 



v(0 = 



v --i 
4 

v 2 



e(£ 



■ti)-—6(t-tD), (11) 
rv 



where <fi(£o) — e and <j>(£i) — e/r. The form of the 
potential for v > v* and small e is sketched in Fig. |J. 
Keep in mind that </> a s(0 is a monotonically increas- 
ing function from e/r at £1 towards the left, and that 
0as(£ - * —00) = 1. As a result, in Fig. |J, V(£) also in- 
creases monotonically towards the left for £ < £1 . On the 
right of £1, V(£) is constant at v 2 /4, and at £o> there is 
an attractive 5-function potential of strength (ru) _1 [ p^[ . 
The crucial feature for the stability analysis below is the 
fact that V(£) stays remarkably flat at a value 2e/r over 
a distance (£1 — £ 2 ) — | lne/r| H, and on the left of £ 2 , 
it increases to the value of w 2 /4 + n — 1, over a distance 
of order unity. 

If there exist negative eigenvalues of the Schrodinger 
equation (|l(i|), then according to Eq. ([)]), r](£,t) grows in 
time in the comoving frame, i.e., the front solution a s(£) 
is unstable. For our purpose, therefore, we look for the 
value of r at which there is a bound state of Eq. ( |l0| ) 
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with eigenvalue E, such that E — > 0— for the poten- 
tial sketched in Fig. ^. This is a problem in elementary 
quantum mechanics. For e — * 0, the potential V(£) is 
essentially constant in the left neighbourhood of £i , and 
hence for v > v* and E — > 0—, ^e(0 can be written as 



= ^ e A 2 («-«o) +Be A 2 («o-C) 



for £ < 6 

for £i < £ < Co 

for£>£o, (12) 



where Ai = ^/w 2 /4— 1 and A2 = u/2. The function 
i/)e{C) must be continuous at £1 and £oj while its slope 
is continuous at £1, but not at £ - Matching of these 
boundary conditions to determine the value of r, where 
the bound state eigenvalue E crosses zero, also requires 
an expression for the distance £0 — £1 . To this end, we 
divide the range of <j> values between and 1 into the three 
regions marked in Fig. |^: (i) region I, where as < s, 
(ii) region II, where e < as < e/r, and (Hi) region III, 
where </> as > e/r. In the comoving frame, the asymptotic 
shape as (£) of the front is the solution of the differential 
equation 4>'^. + v 4>' as + /(</> as ,e,r) = 0, where a prime 
denotes a derivative with respect to £. The solutions of 
0as(£) in the regions I and II that satisfy the continuity 
of 4>as(£) and 4>'as(0i are respectively given by 



<M0 



ee 



-v(t-to) 



TV 



and 
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The length £0 — £1 of region II is obtained by equating 
<Aas(£i) from the second line of Eq. ( |l3| ) to e/r. After 
dividing out a factor of e/r, this condition becomes 



1 



Co -6 



4 = 1. 



(14) 



Thereafter, using Eqs. (|IJ) and (^3|), one arrives at 
Eq. (0). 

The above analysis yields the relation between v and 
the critical value of r in the limit e — > 0. The conver- 
gence with e, i.e., the rate of approach with e of the 
dashed curve to the solid one in Fig. ||, can be obtained 
by considering the effect of the ^</> as _1 (£) term of V(£) on 
the eigenfunctions and eigenvalues. For v > v* , this term 
is simply a correction of order e™ _1 to the finite bottom 
value of the potential. This term can be included pertur- 
batively, and accordingly it leads to a shift of order e n ~ 1 
in the critical value of r. As Fig. || illustrates, this pre- 
diction is confirmed numerically. The case v = v* calls 
for a more detailed analysis, since the bottom value of 
the potential vanishes in the limit e — > 0. In this case, 



that 



,(0~(CC + ^)e-e,Bo V(£) 



it is known [J^ 

n(e/r) n - 1 ((/£x)' n - 1 e( n - 1 ')^-b in the leading order of e. 
In dominant order, we need to keep only the exponential 
behaviour, and the solution of iPe(£,) is then given by 
the Bessel function A^K (2^/ne n - 1 e-( n - 1 ^^-^ /r"- 1 ) 



in the left neighbourhood of £1. The | In -1 e\ scaling for 
the asymptotic approach of the dashed curve to the solid 
one is then easily obtained once the boundary conditions 
at £1 and £0 are matched with the use of Eq. (|l4|). 

The logarithmic convergence of v to v* from below for 
r > r c can be understood from an argument along the 
lines of that for r = 1 Q. For v < v* , the front profile 
0as(£) in region III is of the form <fi as (£,) ~ Csin[fc(£ — 
6) + /?]e~ ? . For r = 1, region II is absent; in that case, 
the matching to the profile in region I and the divergence 
of the width £0 — £2 — | In. ej implies k ~ 7r|ln.£| _1 . For 
r c < r < 1, the matching to region II will change the 
prefactor, but k will still scale as |lne| _1 because the 
width of region III still diverges logarithmically. As for 
r c < r < 1, this translates into a scaling of v* — v as 
I lne|~ 2 , with a prefactor that depends on r. Note that 
this scaling is nicely consistent with the convergence of 
the r-v curve towards the point (r c ,v*) from the left, 
due to the fact that the slope of this curve vanishes at 
this point, and the convergence from below to this point 
scales as the square of the convergence from the left. 

We finally end this paper with the note that if the 
(nonnegative) growth rate is bounded from above by 
f((j),e) in the interval <p < £/ r , but is equal to f(<j>,e) 
for 4> > e/r, then as e — > 0, the asymptotic front speed 
converges to v* with the same logarithmic convergence 
of Eq. (|^) for any r. It simply follows from the inequality 
v as (e/r) < v < v as (e), where v as (e) is given by Eq. (§). 
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This corresponds to 7 = 1 and iV = 10 4 for the model 
considered in 

The 5-function in Eq. (jll|) appears from the functional 
derivative of / in V, since there is a discontinuity of mag- 
nitude e in f(4>asi £> r ) at <j> as — e. This discontinuity con- 
tributes an amount equal to er _1 dO(0 as — e)/d^ as = 
e[r l^asG^ii)!] -1 <K£ - 5o) to V(£)- As as = er according 
to Eq. (|l3|), the prefactor l/(rv). 



